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We consider the transport and the noise characteristics for the case of a T-shape double quan- 
tum dot system using the equation of motion method. Our theoretical results, obtained in an 
approximation equivalent to the Hartree-Fock approximation, account for non-zero on-site Coulomb 
interaction in both the detector and side dots. The existence of a non-zero Coulomb interaction 
implies an additional two resonances in the detector's dot density of states and thereafter affects the 
electronic transport properties of the system. The system's conductance presents two Fano dips as 
function of the energy of the localized electronic level in the side dot. The Fano dips in the system's 
conductance can be observed both for strong (fast detector) and weak coupling (slow detector) be- 
tween the detector dot and the external electrodes. Due to stronger electronic correlations the noise 
characteristics in the case of a slow detector are much higher. This setup may be of interest for the 
practical realization of qubit states in quantum dots systems. 

PACS numbers: 73.63.Kv,72.15.Qm,72.10.-d 



I. INTRODUCTION 

Transport through a complex quantum dot (QD) sys- 
tem is of high interest from both the practical applica- 
tions and theoretical point of view [1[ . From the applica- 
tions point of view, QD systems may provide the perfect 
environment for the implementation of nanoelectronics 
and the realization of quantum bits (qubits). On the 
other hand, QD systems allow the theoretical study of 
quantum many body effects. Particulary, the Anderson 
single impurity model Q was extensively used to under- 
stand the electronic correlations in QD systems. The 
model, in which QD's are represented as impurities, was 
successfully applied to the study of single or multiple-dot 
systems. For example, single QD's systems allow the con- 
trolled realization of the Kondo regime of the Anderson 
impurity problem Q . Multiple-dot systems may be sub- 
ject to inter-dot coupling, and accordingly novel many 
body states can be generated. 

In the case of a double-QD system, with the QD's 
arranged in a series, parallel, or T-shape configuration 
(see Ref. Q for a picture of these configurations), it 
was shown that the arrangement of the component QD's 
plays an important role when transport properties are in- 
vestigated [4( . The results obtained for the system's con- 
ductance can be explained based on Kondo resonances 
influenced by Fano interference effects Q. In the case 
of a series configuration there are no interference effects 
and the possibility of a double peaked Kondo resonance 
due to the inter-dot tunneling may reduce the system's 
conductance. Differently, in the case of a parallel con- 
figuration the two different channels of electron propa- 
gation are responsible for a sharp and a broad Kondo 
resonance. The interference between these two Kondo 
resonances will significantly reduce the system's conduc- 
tance [J]. The double dot T-shape configuration, with 
one dot (detector dot) directly connected to the external 
leads and the second dot (side dot) coupled to the first 



one but not to the external leads, has a particular trans- 
port behavior. In this case, the density of states (DOS) 
of the detector dot (coupled to the external leads) has a 
broad resonance and develops a sharp dip structure due 
to the interference with the states from the second dot 
while the second dot DOS presents a sharp Kondo res- 
onance. Accordingly, the system's conductance is very 
small in the Kondo regime. A similar double dot system, 
with one Kondo dot and one effectively noninteracting 
dot, was shown by Dias da Silva et al. [6[ that it can be 
continuously tuned to create a pseudogapped DOS and 
access a quantum-critical point separating Kondo and 
non-Kondo phases. Different approaches for the study 
of the T-shape system were based on the Anderson im- 
purity model. For example, Wu et al. [7[ considered an 
infinite on-site Coulomb interaction in the detector dot 
so the double occupancy was forbidden in this dot. On 
the other hand, Guclu et al. [8[ assumed an infinite on- 
site Coulomb interaction in the side QD, neglecting it 
in the detector QD. Although the two mentioned con- 
figurations are different, both the Kondo effect and the 
Fano interference effect play an important role in the sys- 
tem's transport properties 0, Q- Transport properties 
in quantum dot systems were discussed also in terms of 
anti-resonance scattering, a method which gives similar 
results with the standard Green's function method 0. 
In short, the evaluation of the system's transmission co- 
efficient shows that for a single impurity Anderson model 
every time the conduction electron's energy is equal to 
the localized level energy, an anti-resonance scattering 
occurs, leading to a dip in the system's conductance. In 
this case the transmission coefficient vanishes and the 
conduction channel is completely blocked by the pres- 
ence of the localized (impurity) level. 

Here we propose an investigation of the T-shape double 
quantum dot system (See Fig. [T]) based on the equation 
of motion (EOM) method. Our approximation is similar 
to the one used by Hewson [l(| in the study of the Ander- 
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FIG. 1: Schematic representation of the T-shape double quan- 
tum dot system. The detector dot (characteristic energy Ed) 
is coupled both to the side quantum dot (characteristic energy 
E a ) and the external electrodes L and R. 



son's single impurity model. We will consider the general 
case with finite on-site Coulomb interaction in both the 
detector and side dots. We will discuss the main trans- 
port properties of the system including the system's con- 
ductance and current noise characteristics. The paper is 
organized as follows: In Section II we present the general 
Hamiltonian of the system and using the EOM method 
we estimate the detector's dot DOS. In Section III we an- 
alyze the system's main electronic transport properties. 
Finally, Section IV presents our conclusions. 



II. MODEL 

The T-shape double quantum dot system is described 
by the following general Hamiltonian: 

H = ^ e k J ktJ . a c ka . a + ^2 E ddld a + U d n d] n di 

k.cr-.a. a 

+ ^2 E a a\a G + U a n a] n ai + t ^ ( d t a v + ald a ) 
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The first term in the Hamiltonian describes the free elec- 
trons in the leads, ct and Cfe CT;a being fermionic cre- 
ation and annihilation operators for electrons with mo- 
mentum k and spin a in the lead a (a = left (L), right 
(R)). The following four terms describe the mesoscopic 
part of the Hamiltonian and correspond to the electrons 
localized in the detector (Ed) and side (E a ) quantum dots 
of the system. Additionally, electrons in each component 
dot are subject to on-site Coulomb interaction described 
by the interaction terms U d and U a , respectively. The 



last two terms in the Hamiltonian describe interactions 
between the system's electrons. The coupling constant t 
characterizes the interaction between the electrons local- 
ized in the detector and side dots of the system. Vkd-a 
characterizes the interaction between the free electrons in 
the lead a and the localized electrons in the detector dot. 
For simplicity we will consider the case Vkd-.L — Vkd;R in 
which the detector couples to the leads only in the sym- 
metric combination Cfc CT = (cko-,L+Ckcr;R)/V2 and the dot 
connects effectively to a single lead, with Vkd = V%VkdL- 
The system's transport properties can be investigated 
using the Green's function formalism. The main quan- 
tity will be the Green's function corresponding to the 
localized electrons in the detector dot. One way to ex- 
tract this Green's function is to use the EOM method. 
It is well known that in the case of a general Anderson 
impurity model, the EOM method leads to an infinite hi- 
erarchy of higher-order Green's functions, so in order to 
obtain the detector's dot electronic Green's function one 
needs to introduce a reliable approximation to truncate 
this hierarchy. The difficulty is mainly introduced by the 
interaction terms in the system's Hamiltonian. When the 
on-site Coulomb interaction term is absent, an exact so- 
lution of the problem is possible as it is well known that 
in this case the set of equations obtained from the EOM 
method are closed. In the case of two fermionic operators 
A and B the Fourier transform of the Green's function 
with respect to the time, Gab(u) = ((A;B)), is given by 
the general equation 



uj((A;B)) = ({A,B}) + (([A,H};B)) , 



(2) 



where (A) represents the mean value of the operator A, 
{A, B} the anti-commutator of the operators A and B, 
and [A, B] their commutator. The last term on the right 
hand side of the equation is responsible for the genera- 
tion of the infinite chain of higher order Green's func- 
tions. Our main goal will be to calculate the (i-electrons 
Green's function, G d d(uj) = ((d a ; dj.)), to be used in the 
estimation of the system's transport properties. The cal- 
culations are relatively simple, however, they lead to a 
chain of coupled Green's functions equations. The ap- 
proximation we used to close this set of equations is sim- 
ilar to the one introduced by Hewson [10( and in some 
sense is equivalent to the standard Hartree-Fock approxi- 
mation. The resulting Green's functions for the localized 
electrons in the detector and side dots are given by 
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where < nn- a > and < n a _<j > represent the average 
occupancy of the electronic levels in the two dots of the 
system. Note that in the case t — the general case of a 
single Anderson impurity is recovered. The above equa- 
tions are coupled as based on the general Green's function 
formalism the average occupancy of the electronic level 
is given by 



{n da ) = 



1 



(5) 



with r] —> 0. A similar relation stands for < n aa >. The 
imaginary part of the Green's functions can be extracted 
if we use the general relation 
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where V represents the principal part and we introduced 
for convenience the notation A = | 1 2 <5(ccJ — £&) 

with 5(x) the delta-Dirac function. By definition the de- 



tector dot DOS is given by 
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(6) 



The simplest situation occurs when the on-site Coulomb 
interaction in both the detector and side dots is neglected 
Ud = U a = 0. In this case, the Green's function are not 
coupled and the detector dot density of state becomes 
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It is convenient to introduce dimcnsionless quantities 
and measure all energies in units of the coupling con- 
stant t relative to the Fermi level of the conduction elec- 
trons in the leads (Ep = 0). We will analyze two differ- 
ent situations relative to the transport properties of the 
T-shape double quantum dot system. In the first case 
A/t > 1 and the conduction electrons will flow relatively 
fast through the detector dot and their interactions with 
the localized electrons in the side dot will not influence 
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FIG. 2: The detector's dot DOS in the absence of on-site 
Coulomb interaction (Ud = U a = 0) for the case of a (a) 
fast detector (A/t = 2) and a (b) slow detector (A/t = 0.4) 
for different values of the ratio E a /t (E a /t — 1.4 - full line, 
E a /t = —1.4 - dashed line, and E a /t = -2.4 - dotted line). 
For both cases E d /t — — 0.2. 



FIG. 3: The detector's dot DOS in the presence of infinite on- 
site Coulomb interaction (Ud = U a = oo) for the case of a (a) 
fast detector (A/t = 2) and a (b) slow detector (A/t = 0.4) 
for different values of the ratio E a /t (E a /t — 1.4 - full line, 
Ea/t = -1.4 - dashed line, and E a /t = -2.4 - dotted line). 
For both cases E d /t — —0.2. 
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FIG. 4: The detector's dot DOS for finite on-site Coulomb 
interaction ((a) and (b) - Ud/t = 0.2, U a /t = 5; (c) and (d) - 
Ud/t = 5, U a /t = 0.2) for the case of a fast detector ((a) and 
(c) - A/t = 2) and for a slow detector ((b) and (d) - A/t = 
0.4) for different values of the ratio E a /t (E a /t = 1.4 - full 
line, E a /t — —2.4 - dashed line). For both cases Ed/t = —0.2. 



Ud = U a = case and use them as a starting point to 
calculate the occupancy of the side dot level < n a a >. 
Thereafter, we use this value to start the self-consistent 
calculation of the two Green's functions using Eqs. ([3]) 
and ([4]). The results of the calculation are presented in 
Figure [3] In the case of infinite on-site Coulomb interac- 
tion the detector dot DOS has the same number of peaks 
as in the previous case, however, in this situation the 
peaks are grouped together. 

A more interesting case is the one with arbitrary on- 
site Coulomb interaction in the two component dots. A 
non-zero on-site Coulomb interaction is responsible for 
two additional peaks in the detector dot DOS. The posi- 
tion of these peaks is determined by the relative strength 
of the on-site Coulomb interaction. Figures St and St> 
considers the detector dot DOS for Ud/t = 0.2 and 
U a /t = 5 and Figures St and SH for U d /t = 0.2 and 
U a /t = 5. Each case is analyzed for E d /t = —0.2 and 
E a /t = 1.4 (full line) and E a /t = -2.4 (dashed line). 
The presence of additional peaks in the detector's dot 
density of states will reflect in the transport properties 
of the system. 

III. TRANSPORT PROPERTIES 



drastically the transport properties. On the opposite sit- 
uation, when A/t < 1, the flow of the conduction elec- 
trons will depend on the interaction with the localized 
electrons in the side dot and therefore we expect this situ- 
ation to reflect on the transport properties of the system. 
We will address the case of a large ratio (A/t = 2) as a 
fast detector and the case of a small ratio (A/t = 0.4) as 
a slow detector. 

In Figure we present the DOS for the detector dot 
for both the fast and slow detector configurations. As 
a general result, the ratio A/t determines the full width 
half maximum (FWHM) value of the density of states. 
One can clearly see that the density of states presents 
sharper peaks for the case of a slow detector. The value 
of the energy level in the side dot (E a ) controls the po- 
sition of the two peaks in the detector's dot DOS. Our 
calculation is performed at T — K , however, we expect 
that effects related to temperature to be minimal on the 
DOS function. For both cases E d /t — —0.2. 

A similar situation occurs for infinite on-site Coulomb 
interaction (Ud = U a = oo). In this case, due to strong 
Coulomb repulsion between electrons, only one electron 
can occupy any of the two localized energy levels in the 
detector and side dots. The two equations giving the elec- 
tron Green's functions in the system's component dots 
are coupled via the terms involving the occupation num- 
bers < nZ > and < n a a >. Although a solution for the 
two Green's functions requires a self-consistent calcula- 
tion, the advantage of the infinite on-site Coulomb inter- 
action case is that the limit U d = U a = oo the form of 
the two equations is simplified. To solve these equations, 
we start from the Green's functions corresponding to the 



Transport properties for the T-shape double quantum 
dot system can be discussed in terms of system's con- 
ductance and current noise characteristics. To calculate 
the system's conductance we use the Meir and Weingreen 
[ll| formula: 



where Go = 2e 2 /h (e is the electron charge and fi the 
Plank constant), pd d {w) is given by Eq.©, and f(x) rep- 
resent the Fermi-Dirac distribution function. This gen- 
eral equation will lead to an expression for the system's 
conductance as function of temperature. We will limit 
our study to the T = K case, when the derivative of 
the Fermi-Dirac function can be substituted by a Dirac- 
delta function. In the case t — 0, i.e., when the side dot 
is decoupled from the system, the general transport the- 
ory for the Anderson single impurity model predicts the 
existence of a Kondo peak in the system's conductance. 
The situation completely changes for t ^ 0. 

The trivial case without on-site Coulomb interaction 
(Ud = U a = 0) can be solved analytically. In this case 
the localized electron's Green's function can be obtained 
exactly and the system's conductance is given by 

A 2 F 2 
G = G - % . (9) 

2tt ( Ed E a - t 2 ) 2 + A 2 £ 2 

When the two dots are decoupled from the leads (A = 0) 
the system's conductance vanishes and when the side dot 
is decoupled from the detector dot (t = 0) the system's 
conductance matches the conductance obtained from the 
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FIG. 5: The system's conductance as function of the energy 
level in the side dot for (a) Ud = U a — and (b) Ud = U a = oo 
in the case of a slow detector A/t = 0.4 (full line) and in the 
case of a fast detector A/t — 2 (dashed line). For both cases 
Edit = -0.2. 



FIG. 6: The system's conductance as function of the energy 
level in the side dot for (a) Ud/t — 5 and U a /t = 0.2 and 
(b) Ud/t = 0.2 and U a /t = 5 in the case of a slow detector 
A/t = 0.4 (full line) and in the case of a fast detector A/t = 2 
(dashed line). For both cases Ed/t = —0.2. 



Anderson single impurity model [12]. The situation is 
slightly different when we consider a nonzero on-site 
Coulomb interaction in both the detector and side dots 
of the system. In this case, an exact analytical solution 
is impossible, approximations being required to obtained 
the localized electrons Green's function. To account for 
the system's transport properties we used the approxi- 
mation introduced by Hewson [Io| which is very similar 
to the Hartree-Fock approximation. The disadvantage 
of this approximation is that it cannot account for the 
physics related to the Kondo effect. Figure [5] presents 
the system's conductance as function of the side dot en- 
ergy level for a) Ud — U a — and b) Ud = U a = oo. 
In both situations when the side dot energy E a is closed 
to the Fermi level (Ep — 0) an anti-resonance scattering 
occurs leading to a sharp drop in the system's conduc- 
tance. As a general result, the dip is sharper in the case 
of infinite on-site Coulomb interaction and a fast detector 
configuration. 

A completely different situation occurs when finite on- 
site Coulomb interaction is considered in both the de- 
tector and side dots. In this case, the system's conduc- 
tance presents a double dip structure corresponding to 
two possible anti-scattering processes corresponding to 
each transport channel. The presence of the two dips 
is consistent and does not depend on the relative size 
of the on-site Coulomb interaction in the two compo- 
nent dots. Figure [5] presents the conduction of the dou- 
ble dot T-shape system for different values of the side 
dot energy level. We consider two different situation, 
Ud/t = 5 and U a /t = 0.2 for Figure [Si and U d /t = 0.2 
and U a /t — 5 for Figure [U In both cases we consider the 
ratio Ed/t = —0.2, however, different values for this ratio 



do not change qualitatively the behavior of the system's 
conductance. 

Quantum and thermal fluctuations are very important 
as they are one of the reasons why quantum correlations 
are difficult to observe experimentally. Such fluctuations 
can be estimated based on the current noise character- 
istics. We will discuss the system's current noise char- 
acteristics in terms of the Fano factor defined as the ra- 
tio of the shot noise S(V) and the current I(V) passing 
through the system, i.e., 7 = S(V)/(2eI(V)), e being the 
electron charge. The current shot noise, S(V), when an 
external bias V is applied to the detector dot is defined 
as a correlation function of current fluctuations and it 
can be proved that is related to the detector dot DOS 
via the transmission function, T(uj) — irApd(u)) [3]: 

4e 2 f eV/2 

S{V) = — <L>T(u,)[l-T(u)], (10) 

n J-eV/2 

where h is the Planck's constant. On the other hand, the 
source current can be calculated as: 

II = T J dco T{u) [f L (u) - f R (cj)} , (11) 

where f a (to) = [exp {(uj- fi a )/k B T} + l] with T be- 
ing the temperature and ks the Boltzmann constant. For 
a symmetrical bias condition we set /i£ = Ep — eV/2 and 
fiR = Ep + eV/2 with Ep — 0. Note that in our partic- 
ular situation the current is conserved, II = Ir- In the 
T = K limit the Fermi-Dirac functions f a (u)) are given 
by usual step functions and the evaluation of the current 
is relatively simple. Special attention should be given to 
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FIG. 7: The system's Fano factor 7 as function of the applied 
bias for (a)-(b) Ud = U a = and (c)-(d) Ud = U a = co in 
the case of a fast detector A/t — 0.4 ((a) and (c)) and in 
the case of a slow detector A/t = 2 ((b) and (d)) for various 
values of the ration E a /t (E a /t = 1.4 - full line, E a /t = — 1.4 - 
dashed line, and E a /t — —2.4 - dotted line). For all situations 
Edit = -0.2. 
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FIG. 8: The system's Fano factor 7 as function of the applied 
bias for (a)-(b) U d /t = 5 and U a /t = 0.2 and (c)-(d) U d /t = 
0.2 and U a /t — 5 in the case of a fast detector A/t = 0.4 ((a) 
and (c)) and in the case of a slow detector A/t — 2 ((b) and 
(d)) for various values of the ration E a /t (E a /t = 1.4 - full 
line, Ea/t = —1.4 - dashed line, and E a /t = —2.4 - dotted 
line). For all situations Ed/t — —0.2. 



the point V — when the Fano factor can be calculated 
as 7 = I-T(Ef). 

In Figure [7] we plotted the Fano factor as function of 
the applied external bias for (a)-(b) = U a = and (c)- 
(d) Ud = U a = 00. We considered both the slow and fast 
detector situations for a fixed energy level in the detec- 
tor dot (Ed/t = —0.2) and different values for the energy 
level in the side dot (E a /t = 1.4 - full line, E a /t = -1.4 - 
dashed line, and E a /t — —2.4 - dotted line). As a general 
result, the Fano factor for the case of a slow detector is 
relatively larger than that of a fast detector. The reason 
for this behavior is due to the difference in the current 
magnitude between the slow and fast detector configura- 
tions, as for both cases the shot noise has the same order 
of magnitude. A similar situation occurs also when we 
consider finite on site interaction in both the detector and 
side dots (see Fig[8]). From the microscopic point of view, 
the larger Fano factor in the slow detector is explained 
by the stronger coupling of the conduction electrons to 
the energy level in the side dot. In all situations the Fano 
factor is relatively larger when the on-site Coulomb in- 
teraction is non-zero, indicating that stronger electronic 
interactions in the system are responsible for more noise. 



IV. CONCLUSIONS 

In conclusion, we presented a theoretical investigation 
of the main transport properties of a T-shape double dot 
system based on the EOM method. Our results are ob- 
tained in an approximation equivalent to the Hartree- 



Fock approximation for zero, finite, and infinite on-site 
Coulomb interaction in both the detector and side dots. 
We considered two possible configurations of the system 
based on the relative interaction strength between the 
conduction electrons and the localized electrons in the 
detector dot. In the slow detector situation, this interac- 
tion is smaller than the interaction between the electrons 
in the detector and side dots (A/t < 1) and we proved 
that the transport properties of the system are strongly 
affected by the presence of the side dot. On the other 
hand, in the case of a fast detector, when the interac- 
tion between the electrons in the detector and side dots 
is relatively small compared to the interaction between 
the conduction electrons and the electrons localized in 
the detector dot (A/t > 1), the role of the side dot is 
diminished in the system's transport properties. The 
first quantity we considered was the detector dot DOS. 
When the on-site Coulomb interaction is set to zero or 
infinite in both dots, the detector dot DOS shows two 
peaks, slightly sharper for the case of a slow detector. 
The situation is different when the on-site Coulomb in- 
teraction has finite value in both the component dots. 
In this case, additional peaks are present in the detector 
dot DOS. Still, in the case of a slow detector, these peaks 
are sharper. The second quantity we considered was the 
system's conductance. As a general result, the system's 
conductance presents an oscillatory behavior as function 
of the energy level in the system's side dot. For the zero 
and infinite on-site Coulomb interaction, there is only one 
oscillation represented by a dip in the system's conduc- 
tance. When the on-site Coulomb interaction is finite in 
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both dots, the system's conductance presents two dips as 
a signature for the additional peaks in the detector's dot 
DOS. As a general result, the dips in the system's con- 
ductance are sharper for the case of a fast detector. This 
situation can be understood if we consider the derivative 
of the Fermi function in the system's conductance which 
can account for all sharp features in the detector's dot 
DOS and lead to a broader conductance for the case of a 
slow detector. Finally, we analyzed the noise in the sys- 
tem's transport in terms of the Fano factor. Our results 
prove that when electronic correlations are stronger the 
noise is higher, making the signal detection in the slow 
detector case more difficult to measure. The Hartree- 
Fock approximation leads to results in good agreement 
with the results generated by the bosonisation method, 
at least when we discuss the system's main transport 
properties. On the other hand, an approximation which 
includes terms beyond the Hartree-Fock approximation 
has to be considered in connection with other possible 
phenomena in quantum dot systems. One example is the 



Kondo effect whose physics can be investigated only if 
we consider additional higher order Green's functions in 
the equation of motion method [HI • Another example is 
the investigation of the differential capacitance in quan- 
tum dots or molecular systems [ljl [l6[ . A higher order 
approximation is required for the understanding of the 
system's transport properties when a transition between 
the Coulomb blockade regime to the Kondo regime is 
considered [l?} . Especially, it will be of great interest to 
account for such effects in multi-quantum dots systems 
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